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Abstract 

We review the ubiquitous presence in multiparticle production processes of quasi-power law distributions (i.e., distributions fol¬ 
lowing pure power laws for large values of the argument but remaining finite, usually exponential, for small values). Special 
emphasis is placed on the conjecture that this reflects the presence in the produced hadronic systems of some intrinsic fluctua¬ 
tions. If described by parameter q they form, together with the scale parameter T ("temperature”), basis of Tsallis distribution, 
f(X) ~ [1 - (1 - , frequently used to describe the relevant distributions (the X being usually a transverse momentum). 

We discuss the origin of such quasi-power law behavior based on our experience with the description of multiparticle production 
processes. In particular, we discuss Tsallis distribution with complex nonextensivity parameter q and argue that it is needed to 
describe log-oscillations as apparently observed in recent data on large momentum distributions in very high energy p-p collisions. 
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1. Introduction 


Multiparticle production processes, which will serve us as the stage to present and discuss the ubiquitous presence 
of the quasi-power law distributions, cover at present a ~ 14 orders of magnitude span in the observed cross sections 
(when one observes transverse momentum, pj distributions UlSlstl)- It came as a surprise that this type of data can 
be fitted by a one simple quasi-power like formula iSIl: 


H{X) = C ■ 


nXQj 


X-n 


for X —> 0, 
for X — > oo, 


( 1 ) 


which smoothly combines pure power-like behavior in one part of phase space with an exponential in another part. 
Before proceeding further, a few words of explanation are in order. Fig. [T] displays our playground, i.e., the phase 
space for particle produced in high energy collision of, say protons, p + p ^ N particles (mostly n mesons). The 
initial energy is and momenta of colliding particles are pA and ps- One usually works in center-of-mass system 
in which, for pp collisions, pA - -pa and \pa,b\ = P)- The momenta of produced secondaries are decomposed in 
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the longitudinal and perpendicular components (with regards to the collision axis), p - [pl,Pt\\ their energies are 

/ ^ -^\l/2 / 9 9\^/^ 

E - Ip + p \ , where p is mass of the produced secondary and p - \p\. Composition [p + \pt\ \ - Pt is called 

transverse mass. It is also customary to use variable y = In [(£ + pi) f (E - p/,)] (rapidity) in which pi = pj sinhy 
and E - pj-coshy. In what follows, we shall be interested in the so called central rapidity region (with y - 0, i.e., 
with pl ^ 0), and in distributions of pt or pj only (i.e., in Eq. ^ X — pj oi: X - pj). 



Figure 1. Schematic view of details of the phase space in multiparticle production. Cf. text for details.. 


Different parts of phase space depicted schematically in Fig. [T] They are dominated by different collision dy¬ 
namics. It is customary to separate them artificially by some momentum scale, transverse momentum parameter 
po, dividing transverse phase space into a predominantly ’’hard” and predominantly ’’soft” part. They distinguished 
themselves by the type of observed spectra of secondaries produced. In the ’’hard” (scarcely populated) region, with 
Pt > Po. they are regarded as essentially power-like, F (pr) ~ Pj", and are usually associated with the hard scattering 
processes between partons (constituents composing nucleons, quarks and gluons) |@1. In the ’’soft” (densely popu¬ 
lated) region, with pj < po, the dominant distribution is exponential one, F(pT)~e'!kpi-pTIT). It is usually associated 
with the thermodynamical description of the hadronizing system with T playing the role of ’’temperature”, with the 
fragmentation of a flux tube with a transverse dimension, or with the production of particles by the Schwinger mecha¬ 
nism 10. However, both formulas can be unihed in a single quasi-power like formula, Eq. ([T]i, interpolating smoothly 
between both regions. It becomes power-like for high pr and exponential-Iike for low pr, as required. One way of 
introducing it is to start from the very large values of pr where we have a pure, scale free, power law. Decreasing now 
Pt towards the demarcation value po and below it, the natural thing to avoid problems with unphysical singularity for 
Pr ^ 0 is to add to pr a constant term, and choose it equal to po. In this way one gets Eq. ([T]i, which for small values 
of X - Pt (and for Xq — po) becomes exponential (Boltzmann-Gibbs - BG) distribution with temperature T = po Q 
This formula coincides with the so called Tsallis nonextensive distribution 0 for n - IKq — 1), 


h,(X) = Cg ■ 






Cq ■ exp^ 



q^l 


Cl • exp 



( 2 ) 


It has been widely used in many other branches of physics 0. For our purposes, both formulas are equivalent with 
n = l/(q - 1) and Xq = riT, and we shall use them interchangeably. Because Eq. © describes nonextensive systems 


^ Formula m is known as the QCD-based Hagedorn formula 0. It was used for the first time in the analysis of UAl experimental data (1 and 
it became one of the standard phenomenological formulas for pj data analysis. 
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in statistical mechanics, the parameter q is usually called the nonextensivity parameter. Eq. (2) becomes the usual 
Boltzmann-Gibbs exponential distribution for ^ ^ 1, with T becoming the temperature. Both Eqs. (1) and (2) have 
been widely used in data analysis (cite, for example, 
processes of multiparticle production (cf., for example 


]J|) and in the phenomenological analysis of 

USEEIMIIIIIII). 

In the next Section |2] we present a review of a number of examples of how such distributions occur (based on 
our experience with applications of Tsallis statistics to multiparticle production processes, but they have general 
applicability). In Section ISTI we present some specific generalization of quasi-power distributions which allows to 
account for a log-periodic oscillations in variable considered. This phenomenon was so far known and observed in all 
other situations resulting in power-like distributions. However, recent multiparticle production data iHSi seem to 
indicate that, apparently, this phenomenon also starts to be visible there and deserves attention. Section 0] contains a 
summary. 


2. How to get a Tsallis distribution - some examples 

In what follows we shall concentrate on Tsallis distributions, Eq. obtained from approaches not based on in 
nonextensive thermodynamic^ We start from examples of constrained systems which lead to Tsallis distributions 
with q < 1. To get also q > I one has to allow for some intrinsic fluctuations (or relaxing constraints). This will be 
shown next. 


2.7. All variables fixed 

In statistical physics, the simplest situation considered is when a system is characterized by variables U - energy, 
T- temperature and N - multiplicity. Usually one or two of them are fixed and the rest fluctuates (either according to 
gamma distribution, in case of U or T, or according to Poisson distribution in case of N, which are integers) ll2^ . Only 
in the thermodynamic limit (i.e., for N co) fluctuations take the form of Gaussian distributions usually discussed in 
textbooks. In 12411 we also discussed in detail situations when all three variables fluctuate inducing some correlations 
in the system. 

However, if all variables are fixed we have distributions of the type of 

with . = (3) 

i.e., Tsallis distributions with q < \. Interestingly enough, such a distribution emerges also directly from the calculus 
of probability for a situation known as induced partition lE^ . In short; N — I randomly chosen independent points 
{Ui, ..., t/jv-i} break segment (0, U) into N parts, the length of which is distributed according to Eq. Q. The length 
of the k'^ part corresponds to the value of energy = Uk+\ - Uu (for ordered Ut). One could think of some analogy 
in physics to the case of random breaks of string in 77 - 1 points in the energy space. Notice that the induced partition 
differs from successive sampling from the uniform distribution, Eu e {0,U - E\ - E 2 - ■ ■ ■ - 7i<r-i], which results in 
/(£) = l/£ 0 . 


2.2. Conditional probability 

To the category of constrained systems also belongs an example of conditional probability. Consider a system of 
n independent points with energies {£,=i,...,jv}i each energy distributed according to a Boltzmann distribution g, (£,) 
(i.e., their sum, U — Ei, is then distributed according to a gamma distribution gN(U)): 


gi (Ei) = i exp 
A 


(-!) 


and 


gdU) = 


1 


(!)-(-!)■ 


A{N - 1) 

If the available energy is limited, U = Na - const, the resulting conditional probability 


f{Ei\U = Na) = 


gi(7',)g^_i(iVa-£,) (N-l) 


gN{Na) 


Na 


1-1^ 

N a 


V-2 


2-q 


l-(l-.)fj 


(4) 


(5) 


^These are discussed in dlli; their theoretical justification is presented in 123] . 
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becomes a Tsallis distribution with 


q = 


N-l> 

N-2 


< 1 


and 


, aN 

A —-. 

N-l 


(6) 


2.3. Statistical physics considerations 

Both above results arise more formally from statistical physics considerations of isolated systems with energy 
U - const and with v degrees of freedom (v particles). Choose a single degree of freedom with energy E « U (i.e., 
the remaining, or reservoir, energy is Er = U - E). If this degree of freedom is in a single, well defined, state then the 
number of states of the whole system is Q(t/ - E) and probability that the energy of the chosen degree of freedom is 
E is P{E) oc EliU - E). Expanding (slowly varying) 


lnQ(t/-£') 


CO 


Z 

*=() 


1 5®lnQ 


with 


1 def 51nQ(£'r) 

~ Wr 


around U, and keeping only the two first terms one gets 


In/’(£’) oc In□(£) oc-jSii, or 


P{E) oc exp(-j8£'). 


(7) 


(8) 


i.e., a Boltzmann distribution (for which q - 1). On the other hand, one usually expects 

(E 






(9) 


(where a\p_ are of the order of unity; we put ai - 1 and, to account for diminishing the number of states in the 
reservoir by one state, a 2 = 2) ll27ll . One can than write that 


dE'i 


= i-lfkl 


0 


<k-\ 


{v-2f 


( 10 ) 


Because 


ln(l+^) = £(-l)' 




k=0 


(k + iy 


the full series for probability of choosing energy E can be written as: 


Q(t/-£ 

P(E) oc-= exp 

Q.(U) ^ 


z 

k=0 


i-iy 


1 


k + 1 (v - 2)^ 


i-m 


it+i 


= c 1- 


v-2 


■f3E 


(v-2) 


= ;8(2-^)[l-(l-^)^£]-. (11) 


For q - l-l/(v-2)< 1 this result coincides with the previous results from the induced partition and conditional 
probability. 

2.4. Systems with fluctuating multiplicity N 

So far, we were getting Tsallis distributions with q < 1. To obtain ^ > 1 one has to relax the restrictions imposed 
on the system considered; for example, by allowing for fluctuations of one of the variables U, N or T. Fluctuations 
of T, known as superstatistics, were considered in |28, 2^ Here we shall consider multiplicity N fluctuating 
according to some distribution PiN). In this case, the resulting distribution is 


where 


f(E) = ^/^(£)W, 

AnIE) - (l “ U - ~ const. 


( 12 ) 


(13) 
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is a distribution for fixed N (to simplify notation we changed - 2 in Eq. (I3]l to N). The most characteristic for 
our purposes distributions P(N) are the, respectively, Binomial Distribution, Pbd, Poissonian Distribution, Ppo and 
Negative Binomial Distributions, Pnbd (cf. 03111 1: 


Pbd(N) 

Ppd(N) 

Pnbd(N) 


N\(_M-N)\\ M ) r 1^) 

N\ 

T{N + k) (<N>\^P I <N>y'‘~^ 

r{N + i)r(k) \ k ) k I 


(14) 

(15) 

(16) 


They lead, respectively, to Tsallis distributions with q ranging from q < \ for the Binomial Distribution, Pbdi Eg. lfTTl l. 
via q - \ Boltzmann distribution for Poissonian Distribution, Ppo, Ea. ifTsT l. to g' > 1 for the Negative Binomial 
Distribution, Pnbd, Ea.(fT9ll (in all cases (3 - {N)IU): 


fBD(E) 

fpD(E) 

fNBD(E) 



exp(-f3E), q^l; 



(17) 

(18) 
(19) 


In all three cases the physical meaning of the parameter q is the same: it measures the strength of multiplicity fluctu¬ 
ations. 


Var{N) _ 1 

< N < N >' 


( 20 ) 


Eor BD one has Var(N)l < N >< I, therefore for it ^ < 1. Eor PD Var(N)l < N >- I, i.e., q' = 1 as well. Eor NBD, 
where Var(N)l < N > >1, one has ^ > 1. 

In the case of ^ > 1, i.e., for NDB, fluctuations of multiplicity N can be translated into fluctuations of temperature 


T. This is possible because, as shown in 015113211 . NBD, which can be written also in the following form. 


P(N)^ 


r{N + k) 

[r(iv+i)r(^)] 


■/(l+r)-^-^ 


can be obtained from the Poisson multiplicity distribution, P{N) - ^e ^, by fluctuating mean multiplicity N using 
gamma distributiorU, 

- k 

r(k) <N > 

Identifying fluctuations of mean multiplicity N with fluctuations of temperature T, one can express the above obser¬ 
vation via fluctuations of temperature. Namely, noticing that 


B - —, < N >- U < B > and y - -^—, 

^ U Hr 

one can rewrite the gamma distribution for mean multiplicity, f{N), as a gamma distribution of mean inverse temper¬ 
ature 0, 


(b) - ^ 1 

f kp \ 

k-\ 

1 expj 

f kp \ 

1 - ' 

(9-1 <I3>) 

r 


■exp| 

f 1 18 \ 

<p> r(/t)' 

[<P>} 

{ <P>] 

iq- 

\)<p>T\ 

[it) 

^ q- \ <p>] 


( 21 ) 


^We have two types of average here: X means average value in a given event whereas <X> denotes averages over all events (or ensembles). 
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And this is precise ly the g amma distribution describing temperature fluctuations, derived, used and investigated in 
superstatistics ^ 281. When convoluted with the Boltzmann-Gibbs distribution, with j3 as scale parameter, it 
carries it into Tsallis distribution with parameter 


q - I + 


Var(fi) 

0 ) 


( 22 ) 


replacing previous Eq. (l20l l and now denoting the strength of temperature fluctuations. 


2.5. Preferential attachment 

So far, to get g' > 1 we were demanding the existence in the system some form of intrinsic (i.e., nonstatistical) 
fluctuations. However, the same effect can be obtained if the system exhibits correlations of the preferential attachment 
type, corresponding to the ”rich-get-richer” phenomenon in networks | ^ 3 351, and if the scale parameter depends 
on the variable under consideration. If 

r ^ = To+ (^-l)£, (23) 

then the probability distribution function, /(£), is given by an equation the solution of which is a Tsallis distribution 
(again, with q > 1 ): 


df(E) 

dE 


1 


T'{E) 


f(E) 


f(E) 


q 


To 


l-(l-q)- 

1 0 


(24) 


For T^(E) - To one gets again the usual exponential distribution. This approach was also applied to an analysis of 
multiparticle production processes in 

The ’’preferential attachment” can be also obtained from superstatistics. As shown above, fluctuations of multi¬ 
plicity N are equivalent to the results of application of superstatistics, where the convolution 


/(£) = ^(7’)exp(-|) 

becomes a Tsallis distribution, Eq. (| 2 ]l, for 

g(T)^ 


dT 


1 


nn)T 


(nTo\" ( nTo\ 


Differentiating Eq. (l25T l one obtains 

df{E) _^ 

dE T{E) 


f(E) where T(E) - To H—. 


(25) 


(26) 


(27) 


This is nothing else but the ’’preferential attachment” case, again resulting in a Tsallis distribution, which for T{E) - 
To becomes BG distribution, cf. Eq. (l24lP[ 


2.6. Multiplicative noise 

Tsallis distribution can also be obtained from multiplicative noise 1129 
equation S, 

^ H- yit)p = f(t). 


defined by the following Langevin 


(28) 


Here 7 (f) and ^(f) denote stochastic processes corresponding to, respectively, multiplicative and additive noises. The 
resulting Fokker-Planck equation for distribution function /, 


dt 


d(Kif) 


dp 


dp^ 


(29) 


^This is not the only place where such a form of T(E) appears. For example, in {23\ it occurs in a description of the thermalization of quarks in 
a quark-gluon plasma by a collision process treated within Fokker-Planck dynamics. 
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where 

K\ = (f> - {y)p and K 2 - Var(^) - 2Cov(^, 'y)p + Var(y)p^, 
has a stationary solution /, which satisfies 


d(K2f) 

dp 


-Kif. 


(30) 

(31) 


If there is no correlation between noises and no drift term due to the additive noise, i.e., for Cov{^, y) - - 0 

the solution of this equation is a Tsallis distribution for p^, 


l+(?-l)y 


fip)- 

If we insist on a solution in the form of Eq. O, 

f(p)- 


■.U rr 2Ear(^) 2Var{y) 

With T = ———; q=\ + 


i + A 
nrj 


<r> 


with n - 


<r> 


q-V 


then the condition to be satisfied has the form; 


K2{p) = 


nT + p 


K\{p)- 


dKxip) 


One then gets a Tsallis distribution (l33T l but now 


n — 2 + 


(y) 


or q - I + 


dp 


Variy) 


Var(y) (y) + 2Var(y) 

and T becomes a ^-dependent quantity (reminiscent of effective temperature Teff as introduced by us in Q): 

T{q) ^{2-q) [To + {q- l)7’i] with To = 

<r> 2<r) 


(32) 


(33) 


(34) 


(35) 


(36) 


2.7. From Shannon entropy to Tsallis distribution 

As shown in S, a Tsallis distribution emerges in a natural way from the usual Shannon entropy, S (for some 
probability density /(x)), by means of the usual MaxEnt approach, if only one imposes the right constraint provided 
by some function of x, h(x): 


S = 


/ 


dxf(x) ln[/(x)] with constraint < h(x) >- I dxf(x)h(x) - const. 


/■ 


(37) 


This approach contains the same information as that based on Tsallis entropy. In fact, one can either use Tsallis 
entropy with relatively simple constraints, or the Shannon entropy with rather complicated ones (cf., for example, a 
list of possible distributions one can get in this way ll^ l. One gets: 


f(x) - exp [To + Ah(x )], 


(38) 


with constants To and T calculated from the normalization of f(x) and from the constraint equation. A constraint 

q-l 


< z >= Zo = 


2 — q 


where z = In 


E 

l-{l-q)- 
1 0 


(39) 


results in a Tsallis distribution (remember that f(z)dz = f(E)dE), 


f(z) = — expl- — 


zo 


zo 


/(£) = 


1 


1 + 


Zo E 


(l+zo)7’o\ l+zoTo 


1 +^ 


2 — q 


\-(l-q)- 
1 0 


(40) 
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To obtain Tq, one has to assume the knowledge of {E) (this would be the an only constraint in the case of BG 
distribution but here it is additional condition to be accounted for). 

Although at the moment there is no clear understanding of the physical meaning of the constraint (l39l l (except 
its obvious usefulness in getting Eq. (l40l) ). it seems to be a natural one from the point of view of the multiplicative 
noise approach represented by Eq. (1281) . That is because there is a connection between the kind of noise occurring 
in Eq. (ijsl l and the condition imposed in the MaxEnt approach 1^ . Namely, for processes described by an additive 
noise, dxidt - ^(f), the natural condition is that imposed on the arithmetic mean, < x >= c + {^)t, and it results in 
exponential distributions. For the multiplicative noise, dxjdt - x-y(t), the natural condition is that imposed on the 
geometric mean, < In .r >= c + {j)t, which results in a power law distribution s. It seems, therefore, that condition 
(l39l) combines both possibilities and leads to a quasi-power law Tsallis distribution combining both types of behavior. 

2.8. Tsallis and QCD 

Recent high energy experiments from the Large Hadron Collider at CERN (CMS 111], ATLAS 12] and ALICE jsll ) 
provided distributions of transverse momenta measured in the previously unprecedented range pj < 180 GeV. The 
measured cross section then spans the range of ~ 14 orders of magnitude. In 0 it was shown that all these data 
can be successfully fitted by Tsallis distribution ([T]i or (|2]i, which is since then widely used in this case (cf.. Fig. |2]as 
example). This caused question, how is it possible because this is the usual domain reserved for the purely perturbative 
QCD approach? In 142 it was demonstrated that: 

• Starting from the pure QCD partonic picture od elementary collisions (’’hard” scatterings between quarks and 
gluons of incoming protons proceeding with high momentum transfer) one gets power index n ^ 4 - 4.5. 
However one observes hadrons which are formed from quarks and gluons by means of complicated branching 
and fragmentation processes. It turns out that all these processes can be parameterized in a relatively simple way 
and one can easily reproduce n ^ 7 - 8 as observed in the experiment (depending on the energy of collision). 

• However, in this way one reproduces properly only the power index n (or q) and resulting distribution is of pure 
power-like type, ~ 1/ pj, diverging for pr 0 instead of being exponential there. To get a Tsallis distribution 
one has to make the same phenomenological step as that proposed to obtain Eq. ([1]): to replace pj pn + pr 
in the denominator of l/pr- So far, the only rationale behind this is that, in the QCD approach, large pr partons 
probe small distances (with small cross sections). With diminishing of pj, these distances become larger (and 
cross sections are increasing) and, eventually, they start to be of the order of the nucleon size (actually it happens 
around pj ~ pro)- At this moment the cross section should stop rising, i.e., it should not depend anymore on 
the further decreasing of transverse momentum pr. This can be modelled by introducing the constant term 
as above. Effectively one then has (llpr) " —» [pro • (1 + PtIpto)] "- The scale parameter pjo can then be 
identified with that in Fig. [1] In a Tsallis fit we use one formula for the whole phase space with pro becoming T 
in the exponent for small transverse momenta, and scale parameter for large pr- Usually one uses an exponential 
formula for pr < pn and power for pr > pro and pjo separates the two parts of phase space. 

3. Log-periodic oscillations 

So far, we have presented possible derivations and applications of relatively simple form of Tsallis distribution 
with Fig. 13 as an example of its apparent success in fitting even the most demanding data so far. However, closer 
inspection of Fig. |3 shows that ratio of data/fit, usually used to estimate the quality of fits, is not flat but shows some 
kind of clearly visible oscillations of log-periodic character, cf. Fig. [3] The first conjunction in such a case is that 
parameters used were not chosen in an optimal way. However, it turns out that these oscillations cannot be eliminated 
by any suitable changes of parameters {q, T) or (m, T) in Eq^ (|2]l or ([T]l, respectively. Here we shall concentrate only 
on data from the CMS experiment |ll], data from ATLAS 12 lead to identical conclusions. One also has to realize 
that to really see these oscillations one needs a rather large domain in pj. Therefore, albeit similar effects can also be 
seen at lower energies, they are not so pronounced as here and will not be discussed here. Assuming that this is not 
an experimental artifact one has to admit that it tells us that the Tsallis formula used is too simple. There still remains 
something hidden in data which has, so far, avoided to be disclosed, and which can signal some genuine dynamical 
effect which is worth been investigating in more detail. 
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Figure 2. (Color online) Fit to large pj data for pp collisions at 0.9 and 7 TeV from CMS experiment using distribution qH. Parameters used 
are, respectively, (T = 0.135, m = 8) and {T = 0.145, m = 6.7). 


Rather than look for another distribution we shall keep to Tsallis distribution and attempt to improve it accordingly 
to account for effects of these log-periodic oscillations observed in data. Because we have two parameters here, power 
index n and scale (’’temperature”) T, the natural approach is to modify one of them. Because data, which we shall 
analyze below, are presented at midrapidity (i.e., for y ~ 0 and longitudinal momentum pi ~ 0) and for large 
transverse momenta, pr » p, the energy £ of a produced particle is essentially equal to its transverse momentum 
Pt, which we shall use in what follows. 


3.1. Quasi-power laws with complex power indices 

First notice that log-periodic oscillations are ubiquitous in systems described in general by power distributions 
1 4411 . Usually they suggest existence of some scale-invariant hierarchical hne-structure in the system and indicate its 
possible multifractality 1451] . In the context of nonextensive statistical mechanics log-periodic oscillations have been 
hrst observed while analyzing the convergence to the critical attractor of dissipative maps and restricted random 
walks II47II . In the case of pure power like distributions the only parameter to manipulated was the power index. It 
was then natural to modify this parameter by allowing it to be complex \ 4^ 4^ . 

For quasi-power Tsallis distributions this idea was first investigated in detail by us in ll4^ . Here we shall apply it to 
multiparticle data as mentioned before. The complex power index results in effective dressing the original distribution 
by multiplying it by a log-oscillating function, usually taken in the form of; 


R(E) - a-\-b cos [c ln(£’ d) f\. 


(41) 


In 114811 we derived such a factor for a Tsallis distribution. For completeness, we recapitulate the main points of 
this derivation. In general, if one deals with a scale invariant function, 0{x), i.e., if 


0(Ax) — pO{x), 


(42) 


then it must have power law behavior. 


0{x) = Cx^"' 


with 


m — — 


Inp 

In 4 


(43) 
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Figure 3. (Color online) Fit to pj dependence of data/fit ratio for results from Fig. [2 Parameters of function R used ED here are, respectively: 
a = 0.865, c = 2.1 for 0.9 TeV and a = 0.909, c = 1.86 for 7 TeV, whereas for both energies b = 0.166, d = 0.948 and / = -1.462. 


It means that, in general, one can write that {k is an arbitrary integer) 




In/r Ink 
In /I In /I 


(44) 


One must now find whether the Tsallis distribution has a similar property and under what conditions. To this end we 
start from differential df{E)ldE of a Tsallis distribution f(E) with power index n, cf. Eq. (|27]i, and write it for finite 
differences llisll . 

6E = a(nT + E), (45) 

where an < 1 is a new parameter. This leads to the following scale invariant relation, 

g[(l + a)x] = (1 - an)g(x) (46) 


where 



This means then that, in general, one can write Eq. ([T]| in the form; 


Six) = 


ln( 1 - an) 2n 

-h ik -. 

ln(l + a) ln(l + a) 


(47) 


(48) 


The power index in Eq. (l48T l (and therefore also in Eq. ([T]i) becomes a complex number, its imaginary part is signaling 
a hierarchy of scales leading to the log-periodic oscillations. 

If we limit ourselves to A: = 0, one recovers the usual real power law solution and mo corresponds to fully con¬ 
tinuous scale invariance. However, in this case the power law exponent mo still depends on a and increases with it 
roughly as 

niQ + —(n + l)a + — h- 3n - l) a^ -h — + 4n^ - n + l^a^ + _ (49) 

The usual Tsallis distribution is recovered only in the limit a —» 0. 

In general one has 


g(x) = ^ Wi • Re (x = X ^ Wk ■ cos [Im (nik) ln(x)]. 

k=0 k=0 
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(50) 
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Figure 4. (Color online) Left panel: The energy dependence of m = mo deduced from CMS (I||andNA49|4l] data and from a compilation [Wibig] 
fol . Right panel: The energy dependence of parameter a present in mo plotted in left panel. 


This is a general form of a Tsallis distribution for complex values of the nonextensivity parameter q. It consists of 
the usual Tsallis form (albeit with a modihed power exponent) and a dressing factor which has the form of a sum of 
log-oscillating components, numbered by parameter k. Because we do not know a priori the details of dynamics of 
processes under consideration (i.e., we do not known the weights Wk), in what follows we only use k — 0 and k — \ 
terms. We obtain approximately, 




1 


Wo -I- Wi cos 


2n 


In 




(51) 


ln(l -H a) 

In this case one could expect that parameters in general modulating factor R in Eq. (HTt could be identihed as follows; 


a - Wo, b-w\, c -- --, d-nT, f--c-ln(nT). (52) 

ln(l + a) 

Comparison of the ht parameters of the oscillating term R in Eq. (HTI) with Eq. (l48l l clearly shows that the observed 
frequency, here given by the parameter c, is more than an order of magnitude smaller than the expected value equal to 
In/ ln(l + a) for any reasonable value of a. To explain this, notice that in our formalism leading to Eq. (|5T i}only one 
evolution step is assumed, whereas in reality we have a whole hierarchy of k evolutions. This results (cf. 11481] 1 in the 
scale parameter c being k times smaller than in (fSTl l. 


2n 

A- ln( 1 + a) 


(53) 


Experimental data indicate that k ^ 22 (for a ^ 0.15 and c ^ 2 ). 

Erom Ea. dTSl l we see that mo > n. This suggests the following explanation of the difference seen between predic¬ 
tion from theory and the experimental data; the measurements in which log-periodic oscillations appear underestimate 
the true value that follows from the underlying dynamics which leads to the smooth Tsallis distribution. As an exam¬ 
ple consider the mo dependence on a, assuming the initial slope n = 4 (this is the value of n expected from the pure 
QCD considerations for partonic interactions 14211 1. The energy behavior of the power index mo in the Tsallis part is 
shown in the left panel of Eig. 0] whereas the energy dependence of the parameter a contained in mo is shown in the 
right panel of Eig. 0] 


3.2. Quasi-power laws with log-periodic scale parameter T 

The phenomenon of log-periodic oscillations observed in data can also be explained in a different way. We can 
keep the nonextensivity parameter q real (as in the original Tsallis distribution) but allow the scale parameter T to 
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oscillate in a specific way, as displayed in Fig. |5] These oscillations can be fitted by a formula similar to Eq. (HTt . 
with generally energy dependent fit parameters (a, h, c, d, /)); 

T = a + b sin [c (ln(£' + + /j. (54) 


5 ’ 

0 ) 

O 


0.155 
0.150 
0.145 
0.140 
0.135 
0.130 
0.125 
0.120 

0.1 1 10 100 
Pj [GeV] 



Figure 5. (Color online) The T = T (pr) for Eq. (3) for which R = 1. Parameters used are: a = 0.132, b = 0.0035, c = 2.2, d = 2.0, / = -0.5 for 
0.9 TeV and a = 0.143, h = 0.0045, c = 2.0, d = 2.0, / = -0.4 for 7 TeV. 


To get such behavior we start from the well known stochastic equation for the temperature evolution 11491] . which 
in the Langevin formulation (allowing for an energy dependent noise, ^(f, £)) has the form: 


dT 1 

— + -T + ^{t,E)T^Q>, 

dt T 

T is relaxation time which, for a while, we keep constant. For the time dependent E - E(t) it reads: 

dT dE 1 


(55) 


(56) 


In the scenario of preferential attachment known from the growth of networks, cf. Section [23] evolution equations 
as given by Eas. (|2^ and (l27l i can be derived from master equation df(E, fjdt - -f{E, t) for the growth of network 
dEjdt given by OsT 

dE E 

-J--- + T (57) 

dt n 

(n coincides with power index in Eq. ([1])). This will therefore be the equation we shall use in what follow^ With it 
one can write Eq. (l56l l as 

(- + t]^ + -T+fit,E)T (58) 

\n / dE T 


^Notice that in the usually multiplicative noise scenario described by y{t), not discussed here, one has ^ = y{t)E + ^(t). 
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This can be subsequently transformed to 


- + Te 
n 


-ln£ 


^Z_ + ir + ^(f,£)7' = <i) 

fl(ln E) T 


and, after differentiating, to 



dT 

2 

\n 1 d(lnEy 

d(\nE) 



Te 


ln£ 


T 


dT 


dat, E) 

+ ^ ^0. 


(i(ln£’) d{\nE) 

For large E (i.e., neglecting terms oc 1 /£) one obtains the following equation for T: 


1 d^T 
nc/(ln£)2 


-+^{t,E) 


dT ^ ^ d^{t,E) 
c/(ln E) d In E) 


= 0 . 


(59) 


(60) 


( 61 ) 


To proceed further one has to specify the energy dependence of the noise ^(f, E). We assume that it increases loga¬ 
rithmically with energy in the following way. 


^(f,£)=^o(0 + —ln£ 
n 


(62) 


(where n is, again, power index from Eq. ([T]) and w is a new parameter). For this choice of noise Eq. dMT) is just an 
equation for the damped hadronic oscillator and has a solution in the form of a log-periodic oscillation of temperature 
with frequency cu: 

T - C exp I -n ■ 


L + llhFl 

2t 2 


In E > ■ sin(a; In £ + <p). 


(63) 


The phase shift parameter <p depends on the unknown initial conditions and is therefore an additional fitting parameter. 
Averaging the noise fluctuations over time t and taking into account that the noise term cannot on average change the 
temperature (cf. Eq. (l55l l in which {dT jdt) — 0 for O = 0), i.e., that 


i+(f(f,£)> = 0, 

r 


(64) 


we have 


T - a -sin(culn£ + (p). 

n 


(65) 


The amplitude of oscillations, b'In, comes from the assumed behavior of the noise as given in Eq. (l62l l. Notice that 
for large n, the energy dependence of the noise disappears. It means that, because, in general, n decreases with energy 
lld^ll . one can expect only negligible oscillations for lower energies but increase with the energy. 

This should now be compared with the parametrization of T (pr) given by Eq. (l54l i and used to fit data in Fig|5] 
Looking at parameters we can see that only a small amount of T (of the order of h/d ~ 3%) comes from the stochastic 
process with energy dependent noise, whereas the main contribution emerges from the usual energy-independent 
Gaussian white noise. 

The above oscillating T needed to fit the log-periodic oscillations seen in data can be obtained in yet another way. 
So far we were assuming that the noise ^(f, E) has the form of Eq. (l62l l and, at the same time, we were keeping the 
relaxation time t constant. In fact, we could equivalently assume the energy E independent white noise, ^(f, E) = 
but allow for the energy dependent relaxation time. Assuming it in the form. 


T = t(E) = 


WTq 

n + co^lnE’ 


( 66 ) 


results in the following time evolution of the temperature. 


T(t) = (T) + [T{t = 0) - {T)]E-“^"/" exp j . 

It gradually approaches its equilibrium value, (T), and reaches it sooner for higher energies. 

13 


( 67 ) 
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3.2.1. Log-periodic oscillations: summary 

To summarize, we have presented two possible mechanisms which could result in the log-periodic oscillations 
apparently present in data for transverse momentum distributions observed in LHC experiments. In both cases one 
uses a Tsallis formula (either in the form of Eq. ([TJ or Eq. (|2]i), with main parameters m - the scaling exponent (or 
nonextensivity q - \ Ijm) and T - the scale parameter (temperature). 


• In the first approach, our Tsallis distribution is decorated by an oscillating factor. This is done by changing in 
Eq. (IZTl i differentials by finite differences, dE 6E - a{nT ■+ E), where the new parameter a < Ijn regulates 
the smallness of dE. As results, we get for x = 1 H-£/(nr) the scale invariant relation, g[(l-i-a)x] = (1 -an)g(x). 
This, in turn, means that power index m (and also nonextensivity parameter q) becomes a complex number, of 
which the imaginary part describes a hierarchy of scales leading to the log-periodic oscillations. The scale 
parameter T remains unaltered. 


It should be mentioned at this point that complex q inevitably also means complex heat capacity C - Ijil - q) 
(c.f., ll^ 3 511). Such complex (frequency dependent) heat capacities (meaning relaxing temperature) are 


widely known and investigated, see 


In the second approach, it is the other way around, i.e., whereas m - n remains untouched, the scale parameter 
T is now oscillating. Erom Eq. (l65l l one can see that T — T(n - 1 /(I -q), E) and as a function of nonextensivity 
parameter q it continues our previous efforts to introduce an effective temperature into the Tsallis distribution, 
Teff = r(^)|[Ill5l 1^ (but here in a much more general form). The two possible mechanisms resulting in 
such T were outlined: (/) - the energy dependent noise connected with the constant relaxation time, or else, (//) 
- the energy independent white noise, but with energy dependent relaxation time. 


At the present level of investigation, we are not able to indicate which of the two possible mechanisms presented 
here (complex q or oscillating T) and resulting in log-periodic oscillations is preferred. This would demand more 
detailed studies on the possible connections with dynamical pictures. Eor example, as discussed long time ago by 
studying apparently similar effects in some exclusive reactions using the QCD Coulomb phase shift idea 0541] . The 
occurrence of some kind of complex power exponents was noticed there as well, albeit on completely different grounds 
than in our case. A possible link with our present analysis would be very interesting but would demand an involved 
and thorough analysis. 


4. Summary and conclusions 


We present examples of a possible mechanisms resulting in the quasi-power distributions exemplified by Tsallis 
distribution Eqs. ([T]i and (|2]l. Our presentation is limited to approaches not derived from nonextensive thermodynamic 
connections of this distribution ifT^ ^ 23 ]. 

It was shown that statistical physics consideration, as well as ’’induced partition process”, results in a Tsallis 
distribution with q - (N - 3)/(N - 2) < 1, Eq. (O. To get q > 1 one has to allow for fluctuations of the multiplicity N. 
They modify the parameter q which is now q - I + Var(N)l{N)^ - 1/{A^), cf. Eq. (l20l) . The conditional probability 
for BG distribution results again in Eq. ([3]). 

We proposed and discussed two possible mechanisms which would allow quasi-power law distributions to which 
Tsallis distribution belongs, to describe data showing a log-periodic ’’decoration” of simple power law distributions. 
One is a generalization of Tsallis distribution to real power n (it can be regarded as generalization of such well known 
distributions as Snadecor distribution (with n = (v -i- 2)/2 with integer v, for v ^ cx) it becomes an exponential 
distribution), can be extended to complex nonextensivity parameter). The other is introducing a specific, log-periodic 
oscillating scale parameter, effective temperature Teff, generalizing our previous results in this field presented and 
used in OIMS]. 

Different derivations of Tsalis distributions turn out to be, in a sense, equivalent. Eluctuations of multiplicity N are 
equivalent to results fluctuations of T (which is the basis of of superstatistics); on the other hand, from a superstatistics 
formula one can get the preferential attachment, cf. Sections [2.4! and l23] 
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